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We investigate the time dependent Ginzburg-Landau (TDGL) equation for a non 
conserved order parameter on an infinitely ramified (deterministic) fractal lattice 
employing two alternative methods: the auxiliary field approach and a numerical 
method of integration of the equations of evolution. In the first case the domain size 
evolves with time as L(t) ~ t 1 / d ™ ) where d w is the anomalous random walk exponent 
associated with the fractal and differs from the normal value 2, which characterizes 
all Euclidean lattices. Such a power law growth is identical to the one observed 
in the study of the spherical model on the same lattice, but fails to describe the 
asymptotic behavior of the numerical solutions of the TDGL equation for a scalar 
order parameter. In fact, the simulations performed on a two dimensional Sierpin- 
ski Carpet indicate that, after an initial stage dominated by a curvature reduction 
mechanism a la Allen-Cahn, the system enters in a regime where the domain walls 
between competing phases are pinned by lattice defects. The lack of translational in- 
variance determines a rough free energy landscape, the existence of many metastable 
minima and the suppression of the marginally stable modes, which in translationally 
invariant systems lead to power law growth and self similar patterns. On fractal 
structures as the temperature vanishes the evolution is frozen, since only thermally 
activated processes can sustain the growth of pinned domains. 
PACS numbers: 64.60C, 64.60M, 64.60A 



INTRODUCTION 

The study of the relaxation dynamics of a system initially in thermal equilibrium and abruptly 
rendered unstable by a sudden change of a controlling field has recently drawn considerable attention, 
not only because many physical properties may depend on the way a material reaches the equilibrium 
state, but also because it poses intriguing problems to the theory such as broken ergodicity, aging 
and dynamical scaling jjj, 

After the quench, i.e. a rapid lowering of the temperature below the critical point of a phase 
separating system, the initial, disordered state looses stability and the system undergoes a coarsening 
process, during which the domains corresponding to different equilibrium phases compete to grow 
in magnitude. About twenty years ago, Allen and Cahn Q realized that when the order parameter 
is non conserved, the driving force towards equilibrium stems from the tendency of the system to 
reduce the curvature of the domain walls and showed that the typical size of the domains , L, 
increases in time with a power law behavior L(t) ~ t 1 !" 1 . In spite of the fact that the theory of phase 
ordering in homogeneous systems is pretty well understood, only recently the ordering kinetics on 
non-translationally invariant, fractal lattices, has become a subject of investigation. 
The present author in collaboration with Petri Q- [0 considered the role of deterministic fractal 
supports, with finite and infinite order of ramification, on the ordering process employing the so- 
called spherical model, which has the advantage of rendering analytical approaches possible. It was 
found by means of an explicit solution that the spherical model on fractal lattices of finite order of 
ramification, such as the Sierpinski gasket of arbitrary embedding dimension does not display a finite 
temperature phase transition; on the contrary on Sierpinski Carpets, whose order of ramification 
is infinite ||- ||, there exist an order-disorder transition provided that the spectral dimension d s 
exceeds the critical value 2 |7J , in accord with the Mermin- Wagner theorem ||] . 
Interestingly, the study of the spherical model with non conserved order parameter has revealed the 
existence, even on fractals, of a characteristic length scale Lit), which increases in time in a power 
like fashion L(t) ~ t 1 ^ , and of dynamical scaling for the correlation functions. Such a dynamical 
exponent z takes the value d w — ln(d + 3)/ln(2) on Sierpinski gaskets of arbitrary embedding 
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dimension d, and d w ~ 2.10 on the planar Sierpinski Carpet, where d w is the random walk exponent. 
These values of z differ from the Allen-Cahn universal value z = 2, which characterizes the diffusive 
domain growth on standard lattices. We have shown |jj that in order to fully characterize the static 
and dynamical properties of the spherical model on fractal lattices two more quantities are required: 
i) the fractal dimension df and ii) the spectral dimension d s , which for many lattices are related to 
d w via the Alexander-Orbach || relation: 

d 

d s 

Unfortunately the analytically soluble spherical model does not yield predictions which can be 
extrapolated to the physically interesting case of the scalar order parameter. In fact, the lack 
of sharp, well defined interfaces between different phases, renders the spherical model physically 
inadequate to describe the phase separation process of an Ising system. 

The purpose of the present paper is to investigate the phase separation process on fractals with 
infinite ramification order, such as the Sierpinski Carpet family , for a non conserved scalar order 
parameter. We expect the latter to be very sensitive to the presence of inhomogeneities in contrast 
with vector fields. The choice of a Sierpinski Carpet has been suggested by the fact that it represents 
perhaps the simplest example of non stochastic fractal lattice, with infinite ramification order || , 
poll . In section II we introduce the Ginzburg-Landau (GL) model on the fractal lattice, in III we 
consider an auxiliary field approach [Tl| | and discuss its asymptotic behavior. In IV we present 
numerical results of the exact equations of motion and monitor in several ways the growth process. 
In V after stressing the similarities between the auxiliary field method and the spherical model, we 
draw the conclusions. 



I. GINZBURG-LANDAU MODEL ON A FRACTAL LATTICE. 

We shall consider a scalar field <p x whose properties depend on a standard Ginzburg-Landau free 
energy functional and defined at every lattice cell, whose coordinate we represent by x: 

D N N 

H[{^}\ = - j E ^ A ^y - 2 E ^ + f E ^ w 

i,j x x 

where r > and g > are the quadratic and quartic couplings of the GL theory while the first 
term is proportional to the surface energy. We shall focus on the dynamical properties of the GL 
model on the two dimensional deterministic Sierpinski Carpet (SC), of fractal dimension (Hausdorff 
dimension) df = In 8/ In 3. 

In order to construct the SC, one divides a square lattice of L x L cells, with L = 3™, into 3x3 
blocks of equal size and the cells contained in the central block are discarded. Dividing again each 
of the remaining blocks into 3x3 sub-blocks and discarding all the central elements as many times 
as necessary to have the smallest sub-blocks constituted of a single cell, one obtains a structure of 
N cells, where N = L d f . 

We make the assumption that the evolution towards equilibrium of the order parameter, <fi x , at the 
site x is given by the Ginzburg-Landau equation: 

— = - T lMtf + ^ t) (2) 

Here rj x (t) represents a Gaussian white noise with zero average and variance : 

< Vx(t)%(t') >= 2T f TS xy S(t - t 1 ) 

where Tf is the temperature of the final equilibrium state, L is a kinetic coefficient and 5 xy the 
Kronecker symbol. 

By substituting eq. ([j]) into (||) we find that (j> x , at any time after the quench, changes according to 
the equation: 

^1 = T[DA xv ^ y (t) + rUt) - 9<Pl(t)] + Vx(t) (3) 
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We shall write A X y as a difference operator in analogy with the discrete representation of the 
Laplacian on Euclidean lattices. Periodic boundary condition are assumed, unless explicitly stated. 
The operator A is defined as: 

A xy = 1 if x, y are nearest neighbor cells 

A xy = otherwise, 
where 2 < Z x < 4 counts the number of nearest neighbors of the site x. 

At temperature Tf = 0, the free energy has two two equivalent minima 4> — ±^/ r /ff- in the case of 
the spherical model we showed that the critical temperature vanishes when the spectral dimension 
is less than 2, while in the scalar case the critical temperature of a Sierpinski Carpet is finite |l2[ . 
In the following, we shall concentrate on the dynamical properties for deep quenches Tf << T c . 



II. AUXILIARY FIELD APPROACH 

In this section the so called auxiliary field method |[| , [jllj and [jl3| , which has provided new insight 
into the ordering dynamics of translationally invariant systems, will be applied to the non conserved 
order parameter dynamics on the Sierpinski Carpet. 

In this approach one replaces the physical field <f> x by an auxiliary field, m x , which varies in a 
smoother fashion across the interfaces and renders approximations feasible. 

One chooses a non linear transformation from <j) x (t) to a new field m x (t) in such a way that the 
latter obeys to an equation simpler than the original one. If such equation is linear, the statistical 
properties of m x are equivalent to those of free gaussian fields and analytical work can be performed. 
Following the presentation of De Siena and Zannetti JTT| ] , one way of determining the transformation 
is to require that the auxiliary field linearizes the local part of the original equation of evolution 
(i.e. a zero dimensional version of eq. ([?]), obtained by setting D = 0). For convenience, unless 
stated explicitly we shall assume the stiffness constant D and the kinetic coefficient T to be 1. The 
auxiliary field m x is introduced via the mapping: 

^ (m * ) = [i + (JW/ 2 (4) 

In order to obtain the equation of motion for the auxiliary field, m x , we need to consider the non 
local term: 

^2 A X y(f>(m y ) = ^2 <X TO y) ~ Z x (j){m x ) (5) 
v <y> 

where the sum ^2 <y> is restricted to the Z x nearest neighbors of the site x. Assuming that <p x is 
a slowing varying function of m Xl one can expand the field at a nearest neighbor site y in a Taylor 
series: 

4>{m y ) = 4>(m x ) + <j)'(m x ){m y - m x ) + -(f)" (m x )(m y - m x ) 2 + H.O.T. (6) 

having indicated with primes the derivatives of <fi x with respect to m x and neglected higher order 
terms (H.O.T.) in the expansion. Collecting together the terms from the Z x nearest neighbors one 
obtains: 

^ A xy (py ~ 4>'{m x ) ^2 A xy m y + -<j)"{m x ) ^ (m x - m v ) 2 (7) 
y y <j> 

The equation of motion for the auxiliary field m x reads: 

= ^v m y + rm x ~ ^Q( m x) /2 ^xyi 171 * ~ m y) 2 ( 8 ) 
y y 

To obtain eq. (ph we have used the identities: 
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g/rcj)(m x f 



4>'(m x ) 

and 

(f>"(m x ) _ 3(g/r)m x 



(9) 



-Q(m x ) (10) 



4>'(m x ) [l + (g/r)m x ] 
To proceed further, we consider a mean- field like approximation for the last term in eq.(||): 

^ A xy (m x - m y ) 2 ~ — A xy (m x - m y ) 2 = -2— ^ A xy m x m y (11) 



Introducing the abbreviation: 

D (t) = -^-^2A xy m x in y (12) 



N ■ 

xy 



the equation of evolution for the auxiliary field m x {t) reads: 

dm x (t) 



dt 



= ^ A xy m y + rm x - Q(m x )D a (t). (13) 



y 



Upon neglecting the last term in eq. (|T3]), i.e. setting Q(m) = 0, we recover the Kawasaki, Yalabik 
and Gunton (KYG) theory juj. Alternatively, we expand the function Q(m x ) ~ 3(g/r)m x to first 
order in the coupling constant g and write pi] : 

= £ A^m^t) + [r - 3^U (i)K(i). (14) 

The field m x thus evolves according to an equation similar to the one found in the spherical model 
and in the Hartree-like approximation, the non linear term being treated self-consistently. 
In order to obtain the properties of the field m x we consider the eigenvalue problem: 

- A xy v a y = e a v a x (15) 

where v x is the x-th component of the eigenvector associated with the eigenvalue e a of the operator 
A. There is an eigenvalue eo = associated with the uniform mode whose eigenvector has all 
elements equal. The asymptotic behavior of the solution of eq. (|l4|) depends on the distribution 
of the smallest eigenvalues. After expanding the field linear superposition of modes of 

amplitude m a : 

N-l 

m x (t) = m a (t)v«. (16) 

at=0 

one finds that each component evolves independently as: 

m a (t) =m a (0)exp[-e a t + B(t)] (17) 
The quantity B(t) must be calculated self-consistently from the governing equation (|i~4|): 

dB(t) g 

—^-=r-3-<D (t)>o (18) 



where the average < • >o is over the initial conditions of the field m x . Using eq. ( [1 2[) and the 
eigenfunction expansion of m x (t) we compute < Do(t) >o as: 

< D (t) > = t a \m a (0)\ 2 exp \-2e a t + 2B(t)}. (19) 
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For N — > oo a continuum density of states approximation p(e) ~ po^ ds ^ 2 1 is appropriate and one 
can write: 

< D (t) >o= A /" de e/9 (e)e- 2e * +2S(t) (20) 

where yC4 is proportional to the amplitude of the fluctuations of the held m x at the instant t = 0. 
One ends with a closed equation for -B(i): 



dB(t) 



= r-i g -e 2B ^A J deep{e)e- 2 ^ (21) 



ft 

From eq. ( |2l| ) one sees that the quantity -B(t) must behave as 

fl(t)~^±*ln* (22) 

Thus asymptotically the quantity < Doit) > goes to a constant value r 2 /(3g), while the equal-time 
correlation function diverges as: 

1 W 

— 22 < m x (t)m x (t) > ~ t. (23) 

a; 

It is possible, now, to compute the correlation function: 

< m x (t)m y (t) > = A J dep(e)e- 2et+2B ^v x (e)v* y (e) (24) 
A scaling form for the above quantity is (see |l5|): 

< m x (t)m y {t) > ~ tey^i-iR/t 1 ' 11 -)^} (25) 

where R = \r x — r y | is the distance between the sites x and y. The above correlation function yields 
the following average value for the droplet size, L, as a function of time: 

<L 2 {t)>ozt^ (26) 

Knowing the evolution of the held m x (t) it is possible to determine the properties of the original 
held 4> x (t). Recalling that m x has a Gaussian distribution at all times, one writes: 

1 r 1 r m t m v „ m x m v , . . 

P m (m x , t-m y ,t) = — cxp {-- -[—f- + —f- - 2 7 — — -} (27) 

Z m l-Y So [x) bo(y) [So{x)S {y)\ L/2 \ 

where Z m represents a normalization factor: 



Z m = 27rVS (x)So(z/)(l- 7 2 ) (28) 

with 

S (x) =< m 2 x (t) > (29) 

and 

G a (x, y; t) =< m x {t)m y {t) > (30) 
v /A (a;)»So(y) 

The average value of the original held x (t) can be calculated from: 
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< 4> x (t) >= / dm x j dm v P m {m x ,t]m y ,t)(l)(m x ) (32) 
while its correlator G(x,y;t) =< 4> x {t)4> y {t) > is given by: 

G(x,y;t) = l dm x I dm v P m {m x ,t]m y ,t)(j){m x )4>{rny) (33) 



Using the form of the distribution function P m and eq.(|33|) one finds for the correlation function of 
the (j> field the result (see jy|): 

< 4>x{t)(j)y{t) >= J dm x J drn v P m (m x ,t;m v ,t)sign(m x )sign(m v ) = ~(t>l oex sm.~ 1 [y(x,y)} (34) 



Substituting eqs.(B3[),(B3) and (61]) into eq. ( 34 ) one obtains: 

< 0.(t)0 y (t) >= -^dn-^exp-CJe/t 1 /*-)^] (35) 

which predicts a smooth decay of the correlation at long distances and the existence of dynamical 
scaling behavior in the late stage growth. 

To check the above results we have solved the equations of evolution ( |l4] ) numerically on a Sierpinski 
Carpet of size 243 x 243 and monitored the growth by measuring the fluctuation of the homogeneous 
component of the order parameter: 



1 N 

c(o,t) = -C£<t> x (t)f 



which is the zero component of the structure factor and grows in time, as shown in fig. (^), with an 
apparent exponent v = 0.93 in agreement with the spherical model result Q which gives v = d s /2, 
with d s = 1.86. The morphology of the field <j) x obtained inserting the solution m x (t) into the non 
linear mapping eq. (|J) is shown in figs.(^j^) and reveals the existence of large droplets growing in 
time in a fashion similar to that observed on compact supports, as if the fractality affects only the 
mass-to-size ratio of the domains, but does not suppress the diffusive motion of the walls. In this 
sense the auxiliary field method agrees quite well with the behavior of the spherical model. 
In the next section we shall compare these findings with a direct simulation of the Langevin equation. 



III. NUMERICAL RESULTS 



We have investigated numerically the non conserved dynamics on the Sierpinski Carpet starting 
from a disordered state, generated assigning to each cell a random number uniformly distributed 
in the interval [—0.125,0.125] and assumed r = g = 1. In order to integrate numerically eq. (@) 
we adopted an Euler discretization scheme with time step At = 0.01 and Sierpinski Carpets of 
different linear sizes L = 27,81,243 with periodic boundary conditions. We have checked that our 
results do not change appreciably if we decrease further At. We also considered several temperatures 
quenches and runs up to t — 1000. The averages for the various quantities presented below refer to 
50 independent random initial configurations. In order to measure the droplet size we have applied 
the Hoshen-Kopelman algorithm |lq| to label the individual clusters formed by nearest neighbor 
cells characterized by the same sign of the order parameter. Quantitative measures of the droplet 
properties are their masses and the radii of gyration. The mass, M, is defined as the number of cell 
belonging to a droplet, while the radius of gyration is defined as: 

1 M 

z 

The sum is over the lattice sites belonging to a droplet, (x z , y z ) is the position of the lattice sites z 
and (x cm , y C m) the position of its center of mass. To compute the mean value of R 2 (t) as a function 
of time we averaged over all droplets and over all initial configurations. 
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The growth law for the average size of the droplets is reported in fig. (^) and the data are compared 
with those referring to an Euclidean square lattice of the same linear size. We observe that the 
growth is much slower in the first system than in the second. In figs. (p|-|rv|) are shown typical 
snapshots of the system: the domain walls, separating opposite phases, sit on locations where the 
surface energy cost is lower and thermal fluctuations are needed to push the system out of these 
minima. Further evidence of the absence of power-law growth stems from the study of C(0, i), which 
is compared in fig. (19) with the corresponding quantity in the Euclidean case, known to grow as 
t d ' 2 . 

Since in the late stage the characteristic width of the interfaces does not change appreciably with 
time, another independent measure of the domain structure is given by the ratio between the total 
number of sites and the number of sites N p covered by interfaces, the so called inverse perimeter 
density 0]: 

P(t) = N sites /N p 

By assigning a perimeter site every time the absolute value of <f> was less than the value 0.75, one 
obtains an estimate of P which is not too sensitive to the above threshold. While P(t) displays 
power law growth in translationally invariant structures and grows proportionally to L(t), in the 
fractal case saturates at a constant value for low temperature quenches as shown in fig. (p0| ) . 
Another quantity with good self-averaging properties is the surface energy density (fig. (|11|)) , whose 
time dependence we studied and compared with the Euclidean case. For compact supports this 
quantity decays as the inverse of the domain size while on the Carpet it displays slow 

relaxation and eventually freezes. 

Based on these results one is lead to the conclusion that, after the early regime, the average domain 
size and other indicators do not grow with a power law behavior, as found in Euclidean lattices and 
in the spherical model on the same fractal support. At zero temperature one finds a breakdown of 
the self-similar dynamical scaling, in contrast with the compact case. Notice that the crossover from 
a diffusive and curvature driven dynamics to a thermally activated, Arrhenius-like, dynamics is not 
disorder induced as in the case of Ising diluted models |l8fl- During the early stage the domains 
coarsen almost linearly with time because the curvature provides the main driving mechanism, while 
at later times the growth becomes much slower and eventually stops since the height of the barriers 
that the domain walls have to surmount increases with the domain size. The dynamics selects 
configurations, in which regions of opposite magnetization are separated by boundaries formed by 
a large number of voids. Since these configurations are associated with low surface tension the 
interfaces remain trapped. 

Interestingly, the crossover time from a power law growth regime, during which the trapping is not 
effective, to a frozen state is independent of the value of the stiffness constant D. This is understood 
by recalling that during the early stage the domain size grows diffusively as: 

R 2 ~ Dt 2/dm 



The pinning typically occurs when a moving interface encounters voids of size of the order of the 
interfacial width w, which scales as hi ~ \[~D. Recalling that the distance between two voids of 
linear dimension w is proportional to w itself, (due to the deterministic nature of the fractal) the 
typical separation of two domain walls is proportional to w. The asymptotic radius of the droplets 
is: 

Rl ~ w 2 ~ D 

which gives a crossover time r: 

Dr Vd w „ D 



Therefore the interfacial stiffness, D determines the typical size of the droplets, but not the crossover 
time scale r at which the freezing occurs. 

Asymmetric, off-critical quenches lead to immediate equilibration of one of the two phases, and 
fragmented patterns are not observed. On the other hand, when one considers the conserved order 
parameter dynamics on percolation structures pi| even off-critical quenches fail to generate the or- 
dered state during a cooling at reasonable rate and lead to frozen dynamic behavior with interrupted 
growth at low temperatures. 
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Twisted boundary conditions. 

We finally looked at the effect on the phase ordering dynamics of a lattice with twisted boundary 
conditions along the x-direction of the lattice, maintaining the periodicity along the y-direction . As 
t — » oo the compact system orders and forms two large domains of the opposite equilibrium phases, 
separated by rather sharp interfaces. On the contrary, on a fractal the system remains frozen in a 
highly fragmented structure as shown in fig.(|l2|) because of the presence of voids, which reduce the 
role of the curvature. To characterize the interfacial roughness we consider the following correlator: 

g{x) = jfYl h x > y)^' y) - (jf y ^ 2 

V y y y 

where N y represents the width of the system and the two arguments denote the longitudinal and 
transverse coordinate respectively. In the case of the regular lattice, the peak of G(x) narrows as the 
interface becomes straighter and sharper, while in the case of the fractal G(x) remains multi-peaked, 
as the interface pervades all the volume as illustrated in fig. (Jl3|). One can go further by defining a 
characteristic width W of the interface as: 

x ^ x X 

As the time t — > oo W(t) saturates to a much larger value than the corresponding quantity in a 
Euclidean lattice with identical boundary conditions and random initial conditions. 



IV. DISCUSSION AND CONCLUSIONS 



On regular lattices, domain walls on sufficiently large scales can be approximated in the continuum 
limit by smooth curved interfaces, which either grow or shrink, whereas in the fractal case the 
background lattice has a pronounced influence on the phase separation. In a translationally invariant 
system the origin of the power law like behavior can be traced back to the interfacial motion: 
immediately after the temperature quench, the system relaxes rapidly along the directions with the 
steepest grad ients of the free energy and the order parameter sets locally to one its equilibrium values 
cp = ±-\/r/g, leaving the system disordered, due to the presence of multiple domains of opposite 
phases separated by thin walls. This early regime is followed by a slower dynamics during which the 
driving force stems from the movement of the domain walls in order to reduce the interfacial energy 
. In other words, a curved portion of an interface will move and its velocity will be proportional 
to its local curvature, while a planar portion in the absence of thermal fluctuations comes at rest. 
The slow behavior of the late stage reflects the flatness of the energy landscape. One can think of 
the representative point of the system as moving along a direction characterized by a very small 
gradient of the free energy 
Let us consider the free energy Hessian : 

x ^ ti = -A xy - [r<p x - 3g<j) x ]6 X y 

d<p X d<py 

its eigenvalue spectrum gives a measure of the stability of a selected configuration and reveals the 
origin of the slow directions of relaxation. If at a given time the configuration is close to one of 
the global minima of H, i.e.. those characterized by uniform values of the field <j) x = ±y/r/g, the 
spectrum contains only positive eigenvalues and the resulting relaxation process is an exponential 
function of time. On the other hand, a symmetric quench leads almost inevitably to the formation 
of a domain structure, which is associated to a saddle point and not to a minimum in phase space 
p2| |. Thus the Hessian calculated in correspondence of the droplet solution must have a negative 
eigenvalue which corresponds to the expanding mode. The approach to equilibrium consists of the 
passage from the unstable droplet configuration to one of the absolute minima. Since the droplet 
size increases in time the free energy becomes flatter and flatter along the unstable path and the 
instability is reduced. Correspondingly the associated negative eigenvalue vanishes giving rise to a 
relaxation slower than any exponential function ]23|] . 

After this premise, we can see how the absence of translational invariance on a fractal support 
modifies all the features of the spectrum of fluctuations described above. 
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• The bulk fluctuations are modified and the system becomes more sensitive to thermal fluctua- 
tions because the density of state exponent d s becomes smaller ( less than d, in fact). However, 
these modes remain stable, since they correspond to positive eigenvalues of the Hessian. 

• The capillary wave spectrum, associated with a kink solution, will also be modified and the 
zero mode Goldstone mode is suppressed by the following mechanism: let us consider a rigid 
shift of a straight interface parallel to one of the directions of higher symmetry of the lattice. 
In fig. ( |l4|) we show the the multivalley structure of the free energy, obtained by changing the 
location of the wall. A sequence of barriers hinders the normal fluctuations of the wall which 
eventually serve to equilibrate the system in a standard lattice. This rugged landscape will 
determine a slower growth of the domains and instead of a diffusive motion in a homogeneous 
space the interface undergoes a diffusion in a system of barriers of varying heights. 

• The negative eigenvalue associated with the presence of a curved interface remains negative 
only for droplets of small size and becomes positive, when the curvature energy becomes of 
the order of magnitude of the barriers. 

Thus the role of the curvature energy, which provides the only driving force in translationally invari- 
ant systems, is greatly reduced, because the lacunarity tends to pin the interfaces in configurations 
which locally minimize the energy. The residual dynamics can be understood in terms of consecu- 
tive transitions among metastable states, i.e. different basins of the free energy. Finally at very low 
temperatures this mechanism eventually ceases and the dynamics becomes completely frozen. 
Let us also remark that the multivalley structures is washed out as the interfacial width exceeds the 
impurity size. Within the limit of infinite thickness one should be able to recover the results of the 
spherical model. In this model, in fact, the width of the walls increases in time during the ordering 
process; as a consequence one observes scaling. 

By contrast, in the scalar model the interfacial width W becomes asymptotically much shorter both 
than the domain size L(t) and the voids, so that the system fails to show power law growth. 
The typical depth B of the local minima of the free energy can roughly be estimated to be propor- 
tional to the size of the domain L to some positive power ip. The value (df — 1) is an upper estimate 
for the exponent ip, since the interfaces do not sit at random, but join places where the surface cost 
decreases. 

Since the late stage growth proceeds only via thermally activated processes one can roughly estimate 
the typical time r for a droplet of size L to overcome a barrier; this is given by the formula: 

r = cxp(B/k B T f ) ~ exp[(L(t)^/k B T f ] 

Thus the domain size L evolves logarithmically as 

L(t) ~ (fcuT/lnt) 1 /* 

in agrrement with the droplet model [ jl8| and stops completely at zero temperature because of the 
prefactor Tf. 

To summarize, we have studied a simple non random system and found that its late stage behavior 
is not characterized by power laws and is not self similar , as in the case of of periodic substrates. 
We have also found that the treatment of the Ginzburg-Landau model on fractal substrates based 
on the auxiliary field method yields results which are consistent with the solution of the spherical 
model on the same support, but does not reproduce the phenomenon of pinning of domain walls 
which is observed numerically. In fact, the spherical model fails to reproduce a necessary feature of 
phase separating systems, i.e. domain walls of finite thickness. Moreover, contrary to the situation 
of the scalar model, within the spherical model there is no activation energy associated with the 
creation of a kink since the energy gap between the ordered phase and the instanton solutions, 
vanishes in the infinite volume limit. This leads in the spherical model to a phase ordering process 
which is diffusive and non-Arrhenius like. The fractality of the substrate has the minor effect of 
rendering less effective the diffusive growth (because d w > 2). Finally the difference between the 
auxiliary field method and the numerical solution seems to be due to the mean-field assumption eq. 
(O), which underestimates the effects of the inhomogeneous lattice. 
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FIG. 1. Peak of C(0, t) (measured in r/g units), employing the auxiliary field method, plotted versus 
time in arbitrary units (a.u.). The line represents the power law growth t da ^ 2 , with d s = 1.86. 




FIG. 2. Morphology of the domains obtained by solving numerically the equations for the auxiliary field 
method on a lattice 243 x 243. Notice that the droplet structure is not strongly affected by the presence of 
the holes, i.e. the interfaces are not pinned. The snapshot refers to times t = 25 in the same arbitrary units 
as fig.(l) 
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FIG. 3. Same as fig. (2), but at t = 250 (a.u.) 




FIG. 4. Configuration of the auxiliary field at t — 2500 (a.u.). 
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FIG. 5. Average size of the droplets (in lattice units) versus time t (a.u.) calculated from the numerical 
solution of the Ginzburg-Landau equation. The data represent the averages over 50 sets of random initial 
conditions and noiseless dynamics T = 0. The upper curve refers to a periodic system 256 x 256, while the 
lower curve to a Sierpinski Carpet of size 243. 




FIG. 6. Instantaneous configurations at t = 10 (a.u.) of the field obtained solving numerically the 
Langevin equation on a Sierpinski Carpet 243 x 243. 
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FIG. 7. Same as fig. (6), but at t = 100 (a.u.) 




FIG. 8. Configuration of the scalar field at t = 1000 (a.u.). Notice that the walls remain frozen at 
positions where the surface energy a local minimum. 
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FIG. 9. Amplitude of the homogeneous component of the order parameter fluctuation versus time (a.u.) 
in the case of non conserved scalar order parameter dynamics on a Sierpinski Carpet 243 x 243 (diamonds) 
and compact square lattice of the same size (crosses). 



FIG. 10. Dimensionless inverse perimeter density versus time in arbitrary units for compact and fractal 
supports. The symbols are the same as in fig. (9) 



s 

S 0.02 



100 200 300 400 500 600 700 800 



FIG. 11. Surface energy in arbitrary units as a function of time for an Euclidean lattice of size 256 x 256 
and for a Sierpinski Carpet 243 x 243. The symbols as in fig. (9) 
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FIG. 12. Typical configuration obtained in the presence of twisted boundary conditions along the vertical 
direction at t = 1000 (a.u.) for a Sierpinski Carpet 243 x 243. 
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FIG. 13. Correlation G(x) computed for for a standard lattice (lines and points) and for a Sierpinski 
Carpet (continuous line) with antiperiodic boundary conditions along the x-direction at time t — 1000 
(a.u.). 
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FIG. 14. Surface energy landscape, associated with a fiat interface as a function of its location, for three 
different choices of the width of the interface. The continuous curve refers to an interface of width w = 4 
in lattice units, the intermediate curve refers to a curve of width w = 8, while the dotted line refers to a 
diffuse (w — 15) interface. Notice the disappearance of small valleys as the interfacial width increases. The 
energy units are arbitrary and the x-axis is in lattice units. 
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